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Introduction 


This  study  considers  the  basic  problem  of  determining  the  flow,  inside  a  free-surface  potential- 
flow  region,  that  corresponds  to  a  given  flow  at  the  surface  bounding  the  potential-flow  region. 
A  new  fundamental  boundary  integral  representation  is  given.  This  representation  of  linear  free- 
surface  potential  flows,  called  velocity  representation,  explicitly  defines  the  velocity  u  within  the 
potential-flow  region  in  terms  of  the  velocity  u  at  the  boundary  surface.  Specifically,  the  velocity 
representation  defines  the  velocity  inside  a  potential-flow  domain  in  terms  of  source  and  vortex 
distributions  with  strength  equal  to  the  normal  component  u  •  n  and  the  tangential  component 
u  x  n  of  the  velocity  at  the  boundary  surface. 

Thus,  the  velocity  representation  does  not  involve  the  velocity  potential  (f>.  This  property  is 
a  major  difference  between  the  velocity  representation  and  the  classical  boundary  integral  repre¬ 
sentation,  called  potential  representation  hereafter,  which  defines  the  velocity  potential  (p  within  a 
potential-flow  region  in  terms  of  the  potential  (j>  and  its  normal  derivative  dcp/dn  at  the  boundary 
surface.  The  velocity  representation  can  therefore  be  used  to  couple  a  viscous  flow,  for  which  a 
velocity  potential  cannot  be  defined,  and  a  potential  flow  in  a  direct  manner,  i.e.  without  having 
to  solve  an  integral  equation  as  is  necessary  if  the  potential  representation  is  used.  Specifically,  the 
potential  representation  requires  a  two-step  procedure:  an  integral  equation  must  first  be  solved 
to  determine  <fr  at  the  boundary  surface  from  the  known  boundary  value  of  d(p/dn ,  and  4>  and  V</> 
can  subsequently  be  determined  at  interior  points  from  the  boundary  values  of  (f>  and  dcpjdn . 

Another  important  difference  between  the  classical  potential  representation  and  the  velocity 
representation  given  here  is  that  the  velocity  representation  defines  u  in  terms  of  first  derivatives  of 
the  Green  function  G ,  whereas  u  can  only  be  obtained  from  the  potential  representation  using  either 
analytical  differentiation  of  the  potential  representation,  which  involves  second-order  derivatives  of 
G ,  or  numerical  differentiation  of  <p. 

The  new  velocity  representation  and,  for  completeness,  the  classical  potential  representation  are 
given  here  for  generic  free-surface  potential  flows  and  for  four  basic  classes  of  flows  corresponding 
to  the  particular  cases  of  (i)  an  infinitely  rigid  or  soft  free-surface  plane,  (ii)  diffraction-radiation  of 
time-harmonic  water  waves  without  forward  speed,  (iii)  steady  ship  waves,  and  (iv)  time-harmonic 
ship  waves  (diffraction-radiation  with  forward  speed).  These  flow  representations,  given  here  for 
deep  water,  can  be  extended  to  uniform  finite  water  depth,  and  indeed  can  be  extended  to  a  broader 
class  of  problems  involving  dispersive  waves  that  propagate  in  the  presence  of  a  planar  boundary. 
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The  velocity  representation  provides  a  new  mathematical  basis  for  computing  free-surface  po¬ 
tential  flows  about  nonlifting  and  lifting  bodies  using  free-surface  Green  functions.  In  particular, 
this  boundary-integral  representation  can  be  used  together  with  the  Fourier-Kochin  approach  ex¬ 
pounded  in  Noblesse  and  Yang  [1] ,  as  is  indeed  shown  here.  Specifically,  the  velocity  representation 
and  the  Fourier-Kochin  approach  are  shown  to  yield  remarkably  simple  analytical  representations 
of  the  waves  generated  by  an  arbitrary  boundary  velocity  distribution  for  time-harmonic  flows,  with 
and  without  forward  speed,  and  for  steady  flows.  The  importance  of  the  Fourier-Kochin  representa¬ 
tions  of  waves  for  practical  purposes  is  demonstrated,  for  the  case  of  steady  ship  waves,  in  Guillerm 
and  Alessandrini  [2]  and  Yang  et  al.  [3,4] .  Specifically,  the  Fourier-Kochin  representation  of  steady 
ship  waves  is  coupled  with  nearfield  calculations  based  on  the  RANS  and  the  Euler  equations  in 
Guillerm  and  Alessandrini  [2]  and  Yang  et  al.  [3] ,  respectively,  and  is  applied  to  the  design  of  a 
wave  cancellation  multihull  ship  in  Yang  et  al.  [4] . 

Coordinates  (f , r? , £)  and  (x,y  ,z)  are  nondimensional  with  respect  to  a  reference  length  L  that 
characterizes  the  size  of  the  ship  or  offshore  structure.  The  velocity  u  —  (u ,  v ,  w)  is  nondimensional 
with  respect  to  a  characteristic  reference  velocity  U,  e.g.  the  forward  speed  U  of  the  ship,  and  the 
velocity  potential  4>  is  nondimensional  with  respect  to  the  reference  potential  U L .  The  2  axis  is 
vertical  and  points  upward,  and  the  mean  free-surface  plane  is  taken  as  the  plane  2  =  0.  For  steady 
and  time-harmonic  flow  about  a  ship  advancing  with  speed  U  in  calm  water  or  in  waves,  the  x  axis 
is  chosen  along  the  path  of  the  ship  and  points  toward  the  ship  bow,  and  the  velocity  u  corresponds 
to  the  flow  disturbance  due  to  the  ship. 

Consider  a  potential-flow  domain  bounded  by  a  closed  boundary  surface  E  defined  as 

£  =  EWUEFUE°°  with  EHU/=E//UElv 

The  surface  Ew  is  an  arbitrary  control  surface  outside  the  viscous  boundary  layer  that  surrounds 
a  ship  hull  (or  other  body,  e.g.  an  offshore  structure)  at  or  below  the  free  surface.  If  viscous  effects 
are  ignored,  may  be  taken  as  the  wetted  surface  of  the  ship  hull.  The  surface  T,w  represents 
the  outer  edge  of  the  viscous  wake  trailing  the  ship,  or  a  control  surface  outside  the  viscous  wake. 
For  a  ship  equipped  with  lifting  surfaces,  e.g.  a  sailboat,  T,w  includes  the  two  sides  of  every  vortex 
sheet  behind  the  ship  hull.  For  a  multihull  ship,  the  hull+wake  surface  T,HW  consists  of  several 
component  surfaces,  which  correspond  to  the  separate  hull  components  of  the  ship  and  their  wakes. 
The  surface  HHW  is  bounded  upward  by  the  mean  free-surface  plane  2  =  0.  The  surface  EF  stands 
for  the  portion  of  the  mean  free-surface  plane  outside  Y,HW.  The  surface  E00  represents  a  large 
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boundary  surface  that  surrounds  T,H .  and  T00  are  the  intersection  curves  between  the  free 

surface  EF  and  the  surfaces  Y>HW  and  E°°. 

Boundary  integral  representations  are  obtained  in  the  study,  for  both  the  potential  (p  and  the 
velocity  u  =  ,  using  elementary  identities  in  vector  calculus  (divergence  theorem  and  Stokes’ 

theorem)  in  the  3D  flow  region  bounded  by  the  closed  surface  E  and  the  2D  free-surface  region  EF 
bounded  by  the  curve  r^^UT00  .  The  surface  E00  and  the  curve  T00  are  not  explicitly  mentioned 
hereafter  because  they  yield  contributions  that  vanish  in  the  limit  when  the  surface  E°°  is  chosen 
infinitely  far  from  the  surface  T,H . 

The  unit  vector  n  =  ( nx,ny,nz )  normal  to  the  boundary  surface  E  points  inside  the  mean  flow 
domain.  Thus,  the  unit  vector  normal  to  the  mean  free  surface  EF  is  n  =  (0,0,— 1).  The  unit 
vector  t  =  (tx, ty, 0 )  tangent  to  the  boundary  curve  r^^UT00  is  oriented  clockwise  (looking  down) 
along  Thw  and  counterclockwise  along  T00.  The  unit  vector  normal  to  the  boundary  curve  rHW 
points  inside  the  mean  flow  domain,  like  n ,  and  is  given  by  V  =  (—ty, tx,0) . 

Potential  and  velocity  representations  for  elementary  Rankine  singularities 

The  velocity  field  u  =  (u,v,w)  at  a  point  £  =  (£,77,0  inside  the  flow  domain  is  given  by 
u  =  Vc  </>(£ )  where  =  (<%  ,  dv ,  c^) .  The  potential  <f>  can  be  expressed  as 

4>  =  ip-X  (la) 

The  potentials  ip  and  x  respectively  correspond  to  distributions  of  sources  and  dipoles  over  the 
boundary  surface  E  : 

ip(£)  —  [ dA(x)  u(x)  -n(x)  G(£,x)  (lb) 

Jn 

x(£)  =  ^cM(z )  4>(x ) VxG(£ , x)-n(x)  (lc) 

Here  Vx  =  (dx ,  dy ,  dz)  and  the  Green  function  G  is  the  potential  of  a  simple  Rankine  source  in 
an  unbounded  fluid,  i.e. 

47tG=—  1/r  with  r  =  \J  (x— £)2+  (y— 7f)2+  (z~  C)2  (2) 

The  velocity  derived  from  the  potential  representation  (1)  is  given  by 

u  =  V^cj)  =  V$ip  —  —  ux  (3a) 
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The  velocity  components  0  and  ux  =  x  are  given  by 


=  J  dA(x )  u(x )  ■  n{x )  V^G(£ ,  x ) 

)  =jdA(£)  0(f)  [ V*G(£ , x)  -nix)} 


The  classical  representation  of  the  velocity  u  given  by  (3)  involves  the  potential  0  in  (3c) .  The  flow 
representation  (3)  is  then  called  the  potential  representation  hereafter.  The  source  component  u'*’ 
is  defined  in  terms  of  the  first  derivatives  of  the  Green  function  G .  However,  the  dipole  component 
ux  involves  second  derivatives  of  G .  Another  representation  of  the  component  ux  that  is  defined 
in  terms  of  a  vortex  sheet  at  the  boundary  surface  and  only  involves  first  derivatives  of  G  is 


)  =  J  dA(x)  [u(f )  x  n(f)]  xV^G(^*,f) 


The  alternative  representation  (4)  of  the  velocity  component  ux  is  identical  to  (3c)  if 


J^dA  [  0  V,£  ( VXG  •  n )  + V^G  x  (Vx  0  x  n )  ]  =  0 


The  relations  V^G  =  -Vx G  and  VXG  =  0,  which  readily  follow  from  (2) ,  show  that  the  x  ,y  ,z 
components  of  (5)  are  given  by  Vx  ■  n ,  Vy  ■  n ,  Vz  ■  n  where  the  vectors  Vx,  Vy,  Vz  are  defined  as 


{4>Gy)y  +  {4>GZ)Z 
—  {(f>Gy)x 
~{<i>Gz)x 


{<f)Gx)y 

vyA  (<t>Gz)z  +  (0Gx)x  >  yz 
I  —  (0G2)y 


-(0Gx)2 

—  {<f>Gy)z 

(0Gx)x  +  ( 4>Gy)y 


We  have  V  •  =  0 ,  V  •  0/  =  o,V-t^=  0.  The  divergence  theorem  applied  to  the  vector  fields 

Vx,  Vy,  Vz  and  the  closed  surface  E  then  shows  that  the  x  ,y  ,z  components  of  (5)  are  null. 

By  substituting  (3b)  and  (4)  into  (3a)  we  obtain 

*(()  =  [  dA{  x)  ([u(f)-n(f)]VeG(£,f)-[u(f)  xn(f)]  xV*G(|\f))  (7a) 

J  s 

which  yields 

u  u-n  G^  +  (uxn)z  Gy  -  (uxn)y  G^  j 

<  v  >  =  f  dA<  u-nGv  +  {uxn)x  Gc  -  (uxn)z  G*  >  (7b) 

J  S 

w  u-n  G^  +  (uxn)y  G^  —  (uxn)x  Gy 

This  expression  defines  the  velocity  u(£ )  at  a  point  £  inside  a  flow  domain,  bounded  by  a  boundary 
surface  E ,  in  terms  of  source  and  vortex  distributions  with  strength  equal  to  the  normal  and  tan¬ 
gential  velocity  components  u-n  and  ux  ft  of  the  velocity  u(x)  at  E  .  The  flow  representation  (7) , 
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which  defines  u(£ )  in  terms  of  the  boundary  velocity  distribution  u(x ) ,  is  called  velocity  represen¬ 
tation.  This  representation  only  involves  the  velocity  u  and  first  derivatives  of  the  Green  function 
G ,  whereas  the  potential  representation  (3)  involves  the  potential  4>  and  second  derivatives  of  G . 

Velocity  and  potential  representations  for  generic  free-surface  flows 

The  potential  representation  (1)  and  the  related  representation  (3)  for  the  velocity  hold  for  a 
generic  Green  function  G  of  the  form 

4? tG  =  -l/r  +  GH 

Here,  the  function  GH  is  harmonic,  i.e.  satisfies  the  Laplace  equation  V2GH  =  0 ,  in  the  region 
bounded  by  the  surface  S .  Green  functions  associated  with  boundary  conditions  at  the  mean 
free-surface  plane  z  =  0  can  be  expressed  as 

G((,x)  =  Gs(x  -  £  ,y  -  tj ,  z  -  Q  +  GF (x  -  £  ,y  -  rj ,  z  +  Q  (8a) 

Gs  corresponds  to  a  simple  Rankine  source  and  is  given  by 

47rG5=-l/r  with  r  —  yj  (x— £)2+  (2/-t?)2+  (z  — C)2  (8b) 

Gf  accounts  for  free-surface  effects  and  is  harmonic  in  the  region  bounded  by  E .  The  simple- 
singularity  component  Gs  is  a  function  of  the  three  variables  (x  —  £,y  —  rj  ,z  —  £) .  The  free-surface 
component  GF  is  a  function  of  (x  —  £,y  —  rj,z  +  Q  and  satisfies  the  relations 

( Gf,  Gf,Gf)  =  {-GFx,  -Gf,  Gf )  and  VX2GF  -  0  (9) 

The  potential  representation  (3)  and  the  decomposition  (8a)  of  the  Green  function  G  show  that 
the  velocity  field  u  can  be  expressed  as 

u  =  us+uF  (10a) 

The  components  us  and  uF  correspond  to  Gs  and  GF  in  (8a) .  Expressions  (8b)  and  (2)  show 
that  the  simple-singularity  component  us  can  be  expressed  in  the  form  (7b) ,  i.e. 

— u-n  Gf  —  (uxn)z  Gf  +  (• uxn)y  Gf 
-u-n  Gf  —  ( uxn)x  Gf  +  ( uxn)z  Gf 
-u-n  Gf  -  (1 uxn)y  Gf  +  (u x n)x  Gf 
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It  is  shown  in  Appendix  1  that  the  free-surface  component  u  F  can  similarly  be  expressed  as 

,F 


IT 


dAl 


—w‘ 


-u-n  GF  —  (ux n)z  Gy  +  (ux  n)y  Gf 
-u-n  Gy  —  (ux n)x  Gf  +  (ux n)z  Gf 


(10c) 


-u-n  Gf  —  (ux n)y  Gf  +  (ux n)x  Gy  j 
The  representations  (10b) ,  (10c)  can  be  expressed  in  compact  vectorial  forms  similar  to  (7a) 

us  =  [  dA[(uxn)  xVxGs-(u-n)VxGs] 

(uF,  vF ; -wF)  =  f  dA  [  (u  x  n )  x  VxGf  -  (u  ■  n )  VXGF  ] 

Jt, 


(lOd) 

(lOe) 


Expressions  (10) ,  (1) ,  (8a)  yield 


N 

u 

>  =  /  dA< 

V 

w 

k 

u-fi  G+  -<t>  (nx  G+  +  nyG+  +  nz  G+) 
-u-n  Gf  —  (uxn)z  G+  +  (uxn)y  G+ 
-u-fi  Gf  —  (uxn)x  Gf  +  (uxn)z Gf 
[  -u-n  G~  -  (uxn)y  G~  +  (uxn)x  G~  J 


(Ha) 


where  G±  is  defined  as 


(lib) 


G±  =Gs±Gf 

The  free-surface  component  GF  in  (8a)  and  (lib)  can  be  expressed  in  the  alternative  forms 

l/r'  +  H 


4t tGf  = 


-l/r'  +  H 


(12a) 


with 


r'  =  y/  (a:-£)2  +  (y-r?)2  +  (*+<)2 

Expressions  (lib) ,  (8b)  and  (12a)  yield  the  alternative  representations 

r 

,  —l/r  +  l/r'  +  H 
47r  G±  =  < 

-l/r  +  l/r'  ±H 

< 

By  substituting  (13)  into  (11a)  we  can  express  0  and  u  in  the  alternative  forms 


(12b) 


(13) 


1 

+ 

* 

>  U  =  < 

uR  -  uR  +  uH 

^cf)R  +  +  4)n 

y  uR  +  uR  +  un 

(14) 
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The  components  <j>R ,  uR  and  <f>R ,  uR  are  defined  by  distributions  of  elementary  Rankine  singular¬ 
ities  1/r  and  1  /r' .  Specifically,  the  Rankine  component  <f>R ,  uR  is  given  by 


'-<j>R 

-(u-n)R  +  4>  ( nxRx  +  nvRy  +  nzRz) 

47T  < 

UR 

>=  dA< 

-(u-n)Rx  —  ( uxn)z  Ry  +  ( uxn)y  Rz 

VR 

Jt, 

-{u-n)Ry  -  (uxn)x  Rz  +  (uxn)z  Rx 

WR 

\  J 

-(u-n)Rz  —  (uxn)y  Rx  +  ( uxn)x  Ry  . 

with  R  =  —1/r .  The  Rankine  component  <f>R .  uR  is  given  by 


' -(u-n)R*  +  4> ( nxR*x  +  nyR*y  +  nzR*z) ' 

An  < 

uR 

*  =  I  dAi 

-{u  ■  n )  R*x  -  (u x n)z  R*y  +  (ux n)y  R*z 

VR 

Jt. 

-(u-n)R*  -  (uxn)x  R*z  +  (uxn)z  R*x 

—wR 

<  *  J 

-{u-n)R*z-{uxn)y  R*x  +  {uxn)x  R*y  j 

with  R*  =  —1/r' .  The  components  <f>H,  uH  and  (j)R,  un  are  defined  in  terms  of  the  functions 
H  and  “H ,  which  account  for  free-surface  effects.  The  free-surface  component  (j/7 ,  uF  —  where  T 
stands  for  either  H  or  T~L  —  is  given  by 


-{u-n)F  +  <j>  (nxFx  +  nyTy  +  nzTz) 

47 r  < 

w r 

>  =  /  dA  < 

-{u-n)Fx  -  (uxn)z Ty  +  (uxn)y Tz 

VT 

Jt. 

-{u-n)Ty  -  (uxn)x Tz  +  (uxn)z Tx 

—vf7 
\  ^ 

— (u  •  n')Tz  —  (uxn)y  Tx  +  (uxn)x  Ty  * 

The  flow  representation  (14)-(16)  is  valid  for  a  generic  Green  function  G  of  the  type  defined  by 
(8)  and  (12) .  This  representation  defines  a  free-surface  potential  flow  in  terms  of  the  flow  at  a 
boundary  surface  £ ,  which  includes  the  free-surface  plane  HF .  The  contribution  of  T,F  is  now 
considered. 


Free-surface  contribution  for  generic  free-surface  flows 

The  contribution  of  the  mean  free  surface  £F,  where  2  =  0  and  (n1,  ny,  nz )  =  (0 , 0 ,  — 1) ,  to  the 
potential  and  the  velocity  (11a)  is  given  by 


—w  G+  +  (j)  GF 

U 

>  =  /  dA< 

wG+  +  uGf 

> 

V 

Jt;f 

wGF  +  vGf 

w 

w  Gz  —  uG~  —  v  G~ 
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At  the  free-surface  plane  2  =  0 ,  (13)  yields 


4t r  (G+,  G+,  G+  G2+)  =  (tf ,  Hx  ,  ^ ,  W«) 

4tt  (G-  G;,  G-,  G~)  =  -(H ,  Hx  ,  Hy ,  tf2) 

These  relations  show  that  the  contribution  of  the  mean  free  surface  HF  is  given  by 

—w  H  +  <t>Hz 
wHx  +  uHz 
w  Hv  +  vHz 
uHx  +  vHy-wHz 

The  free-surface  contribution  (17)  is  considered  below.  The  two  simple  special  cases  of  an  infinitely 
rigid  or  infinitely  soft  free  surface  are  examined  first.  Three  basic  classes  of  free-surface  flows, 
corresponding  to  diffraction-radiation  of  time-harmonic  water  waves  without  forward  speed,  steady 
ship  waves,  and  time-harmonic  ship  waves  are  considered  subsequently. 


Special  cases  of  rigid  and  soft  free-surface  planes 

The  free-surface  contribution  (17)  is  null  if  the  boundary  condition  at  T,F  is  w  =  0  and  the 
function  H  in  (12a)  is  chosen  as  H  =  0,  so  that  the  Green  function  satisfies  the  condition  G^  =  0 
at  (  =  0 .  The  free-surface  contribution  (17)  is  also  null  if  the  boundary  condition  at  EF  is  4>  =  0 , 
which  implies  u  =  0  =  v,  and  the  function  H  in  (12a)  is  taken  as  H  =  0,  which  yields  G  =  0  at 
£  =  0 .  Thus,  the  potential  <p  and  the  velocity  u  are  defined  by  (14)  as 


^  4>R  +  <I>*  f 

(f)F  —  <f>F 

<  j 


at  z  —  0 


(18) 


Here,  4>R,  uR  and  ,  uR  are  given  by  (15)  where  the  boundary  surface  E  is  taken  as  the  hull+wake 
surface  T,HW . 


Wave  diffraction-radiation  without  forward  speed 

The  free-surface  contribution  (17)  is  now  considered  for  the  free-surface  boundary  condition 

w  =  f‘2(j)  +  p  at  2  =  0  (19) 

associated  with  diffraction-radiation  of  time-harmonic  water  waves  without  forward  speed.  Here, 
f  =  u  \/L/g  is  the  nondimensional  wave  frequency,  and  p  stands  for  a  nonhomogeneous  forcing 
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term.  E.g.,  p  may  account  for  a  pressure  distribution  at  the  free  surface.  The  corresponding  Green 
function  G  —  Gs+  GF  satisfies  the  free-surface  boundary  condition 

Gc  =  /2G  at  C  =  0  (20a) 

The  decompositions  (12a)  then  yield  =  f2H  at  C  =  0 .  This  relation  also  holds  for  (  <  0  since 
H  and  H  are  functions  of  2  -f  C, .  We  then  have 


fH  =  nc  =  nz 

It  follows  that  the  Green  function  also  satisfies  the  condition 

Gz  =  f2G  at  2  =  0 

By  using  (19) ,  (20b) ,  and  the  Laplace  equation  V2H  =  0  in  (17)  we  can  obtain 


(20b) 


(20c) 


f  \ 

<t> 

f 

-pH 

47T  < 

u 

>=  dA< 

{(t>'Hz)x  +pHx 

V 

JzF 

(4>  H^y  +  P  Hy 

w 

K  > 

„  (4>%c)x  +  {<t>Hy)y-pHZi 

The  flow  representation  (14)  then  yields 

(j)  =  (f)R  +  (f>R  +  (jfl  —  (f>P 
u  =  uR  +  uR  +  un  +  uV-  —  up 


(21a) 


The  components  <f>R,  uR  and  <f>R,  uR  in  (21a)  are  given  by  (15)  where  the  boundary  surface  E  is 
restricted  to  the  hull+wake  surface  HHW .  The  component  (jfl,  uR  is  given  by  (16)  with  'Z  —  T,HW 
and  T  =  H .  Stokes’  theorem  shows  that  the  component  u^,  which  accounts  for  the  free-surface 
contribution  (17)  if  p  =  0,  is  given  by 


(21b) 


Finally,  the  component  (j)p.  up  associated  with  the  nonhomogeneous  term  p  in  the  free-surface 
condition  (19)  is  given  by 


(-?1 

ty  Hz 

IvP 

>  =  [  dC  (j)  < 

-txHz 

J  rHW 

KJ 

tyHx-txHy  f 

47T  < 


\ 

V 

f  \ 

H 

H 

up 

>  =  /  dA  p  < 

-Hx 

>=  dA  p  < 

H* 

vp 

Jy,f 

~Hy 

Jt.f 

Hy 

wp 

Hz 

\  / 

H< 

(21c) 
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(8)  and  (12)  yield 


II 

s 

M 

I UP j  J*F 

The  representation  (21a)  consists  of  two  components,  which  define  the  potential  4>  and  the  ve¬ 
locity  u .  The  potential  component  of  (21a)  is  a  classical  representation  that  defines  4>  in  the  flow 
domain  in  terms  of  the  values  of  u  ■  n  and  <j>  at  the  boundary  surface  Y,HW .  This  representation 
involves  both  the  Green  function  and  its  first  derivatives.  The  velocity  component  of  (21a)  is  a  new 
boundary-integral  representation  that  defines  u  in  the  flow  domain  in  terms  of  u  at  the  boundary 
surface  '£HW  and  0  =  (■ w-p)/ f 2  at  the  boundary  curve  THW .  The  representation  of  u  only  involves 
VG ,  i.e.  it  does  not  involve  either  the  Green  function  G  or  higher-than-first  derivatives  of  G .  The 
velocity  component  of  (21a)  involves  the  line  integral  (21b)  around  the  boundary  curve  YHW \  the 
potential  component  of  (21a)  does  not  involve  such  a  line  integral. 

Steady  ship  waves 

The  free-surface  contribution  (17)  is  now  considered  for  the  Kelvin  free-surface  condition 

w  +  F2ux  =  p  at  2  =  0  (22) 

associated  with  steady  ship  waves.  Here,  F  =  U/\fgL  is  the  Proude  number,  and  p  stands  for  a 
nonhomogeneous  forcing  term  as  in  (19) .  The  Kelvin  source  potential  G  =  Gs+  GF  satisfies  the 
free-surface  boundary  condition 

GC  +  F2G^  =  0  at  C  =  0  (23a) 

The  decompositions  (12a)  then  yield  Hq  +  F2H ££  =  0  at  £  =  0 .  This  relation  also  holds  for  £  <  0 

because  H  and  H  are  functions  of  2  +  £ .  We  then  have 

=  - F2H (23b) 
It  follows  that  the  functions  H  and  H  and  the  Green  function  G  also  satisfy  the  relations 

Hz  =  -F2Hxx  H  =  -F2Hzxx  (23c) 

Gz  +  F2Gxx  =  0  at  2  =  0  (23d) 

The  function  Hz  in  (23c)  is  defined  as 

Hz(x-£,y-r),z  +  C)  =  /  dt  H(x  -  £,y -r)  ,t  +  0  (24) 

J  —  OO 
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By  using  the  free-surface  condition  (22) ,  the  relations  (23c) ,  the  irrotationality  condition  vx  =  a 
and  the  Laplace  equation  V2HX  =  0  in  (17)  we  can  obtain 


> 

<t> 

<t>p 

u 

>  =  -F2[  dA  Vf-Af< 

up 

47 r  < 

V 

vp 

w 

K  > 

wp 

where  <f>p  and  up  are  given  by  (21c)  and  Vp  is  the  vector  defined  as 


VF 


{(j)  u  H)x 
(u  Ex)x 

{v  Hx)x  (u  Hx)y  (u  Hy)x 

k  {vHzxy)x  -  (uH*y)v  -  (uHz)x  j 


The  flow  representation  (14)  then  yields 


(25) 


<t>  =  (f)R  -  +  (j)H  +  (j)**  —  (f)P 

U  =  UR  —  UR  +  UH  +  Up  —  up 

The  components  </>R,  uR  and  d>R,  uR  in  (26a)  are  given  by  (15)  where  the  boundary  surface  £  is 
restricted  to  the  hull+wake  surface  Y,HW .  The  component  <pH,  uH  is  given  by  (16)  with  £  =  T,HW 
and  T-E.  Stokes’  theorem  shows  that  the  component  <j>R ,  Up  associated  with  the  free-surface 
contribution  (17)  is  given  by 


'utyH-4>tyHx 

uR 

o  r 

utyHx 

47 r  < 

1 

II 

Jthw 

u tyHy  +  u-t  Hx 

-mH 
Wp  j 

u tyHz  —  u-t  HI 

The  identity  u  =  (u  ■  t )  tx -  (u  -  v)  ty ,  where  t  =  (tx,  ty)  and  u  =  (~ty,  tx)  are  unit  vectors  tangent 
and  normal  to  the  boundary  curve  rHW,  yields  the  alternative  representation 


1 

V. _ 

'  (u  ■  v)  (ty)2H  -  ( u-t)txtyH  +  <j>tyHx  ' 

An  < 

Up 

- 

CN 

II 

(u  ■  v)  (ty)2Hx  -  (u  -t)txtyHx 

VF 

JThw 

(u-u)  ( tyfHy  -  {U-t)  {tXtyHy  +  Hx) 

-wH 

wF  j 

M-v)  C ty)2Hz  -(u-t)  (txtyHz-H*xy) 

Both  the  potential  and  velocity  components  of  (26a)  involve  a  line  integral  around  THW .  The 
potential  component  is  a  classical  representation  that  defines  0  in  the  flow  domain  in  terms  of  the 
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values  of  u  ■  n  and  <p  at  Y,HW  and  the  values  of  u  and  (p  at  .  This  representation  involves  G  and 
VG.  The  velocity  component  of  (26a)  is  a  new  representation  that  defines  u  in  the  flow  domain  in 
terms  of  the  boundary  values  of  the  normal  velocity  u  ■  ft ,  u  ■  u  and  the  tangential  velocity  u  x  n, 
u-t.  Thus,  the  velocity  representation  does  not  involve  <fi .  This  representation  involves  VG  and 
the  function  H£y  defined  by  (24) ,  which  is  comparable  to  a  first  derivative  (H*y  and  Vx H  have 
similar  behaviors  in  the  near  field  and  the  far  field) . 

Time-harmonic  ship  waves 

The  free-surface  contribution  (17)  is  now  considered  for  the  free-surface  condition 

w  -  f2(p  +  i2TU  +  F2ux  =  p  at  z  =  0  (27) 

corresponding  to  wave  diffraction-radiation  with  forward  speed.  Here,  /  and  F  are  the  nondimen- 
sional  wave  frequency  and  the  Froude  number,  and  r  =  fF.  The  ship-motion  Green  function 
G  =  Gs+  Gf  satisfies  the  free-surface  boundary  condition 

Gc-/2G  +  z'2tGc-(-F2G^  =  0  at  C  =  0  (28a) 

The  decompositions  (12a)  yield  H q  —  f2H  +  i  2 r  +  F2H ^  —  0  at  (  =  0 .  This  relation  also  holds 
for  C  <  0  since  H  and  H  are  functions  of  2  +  C  •  We  then  have 


=  f2H  -  i  2 r  Hs  -  F2H £  e  (28b) 

It  follows  that  the  functions  H  and  H  and  the  Green  function  G  also  satisfy  the  relations 

7iz  =  f2H  +  i2tHx-F2Hxx  H  =  f2Hz  +  i2tHzx-F2Hzxx  (28c) 

Gz  -  f2G  -i2rGx  +  F2Gxx  =  0  at  z  =  0  (28d) 

Hz  in  (28c)  is  defined  by  (24) .  By  using  the  free-surface  condition  (27) ,  the  relations  (28c) ,  the 
irrotationality  condition  vx  =  uy  and  the  Laplace  equation  V2/P  =  0  in  (17)  we  can  obtain 


4> 

( 

V 

u 

>  =  [  dA(f2Vf  +  i2TVT-F2VF)~4ir< 

uP 

47T  < 

> 

V 

Jef 

v*> 

w 

\  > 

wp 

\ 
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where  the  component  (j)p,  up  is  defined  by  (21c)  and  the  vectors  VF,Vf,  VT  are  given  by  (25)  and 


N 

0 

(<i>  H)x 

{4>H)X 

>  vr  =  < 

0 

.  w  fffl, +w  , 

The  flow  representation  (14)  yields 


j  (f>  =  (pR  -  (f)R  +  (j)H  +  (f)H  +  (j)H  _  0P 

(  U  =  UR  —  UR  +  UH+  Up  +  uR  +  UR  —  up 


(29a) 


The  components  (j>R,uR  and  4>R,  uR  in  (29a)  axe  given  by  (15)  with  T,  =  EHW.  The  component 
<t>H ,  uH  is  given  by  (16)  with  E  =  Y>HW  and  T  =  H.  The  component  <pR  ,  uR  is  given  by  (26b) . 
Stokes’  theorem  shows  that  the  components  d>R  ,  uR  and  uf  are  given  by 


47T  < 


<ptyH 

UR 

f 

0 

>=i2r  dC  < 

> 

vth 

JrHW 

u  -t  H 

wth 
•  > 

u-tm 

<  y 

(29b) 


vH  ) 
f 

PH 

47 r  < 

V?  \  =  f2  dC  (j)  i 

—txH  > 

J  I  JrHw 

L-/J 

tvm-txm 

<  y  a 

In  the  special  cases  /  =  0  and  F  =  0,  (29)  is  identical  to  (26)  and  (21)  as  one  expects. 


(29c) 


Fourier  component  in  free-surface  Green  functions 


Free-surface  effects  in  the  potential  and  velocity  representations  (21a),  (26a),  (29a)  are  defined 
in  terms  of  the  components  ft  or  H  in  the  alternative  decompositions  (12a)  of  the  free-surface 
component  GF  in  the  representation  (8a)  of  the  Green  functions  associated  with  the  free-surface 
conditions 

f  w  -  f2<p  =  p  } 


w  +  F2ux  =  p 


at  2  =  0 


(30) 


w  —  f2(f)  +  i2ru  +  F2ux  =  p 
for  diffraction-radiation  of  time-harmonic  waves  without  forward  speed,  steady  ship  waves,  and 
diffraction-radiation  of  time-harmonic  waves  with  forward  speed.  The  free-surface  components  ft 
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and  H  are  given  by  the  Fourier  superposition  of  elementary  waves 


with 


T=  { 


F= 

1  r oo  roc 

-  lim  —  dp  c 

~  >+0  7 T  J —oo  J —oo 

-] 

r  } 

n  ] 

[  f/k 

j 

>  A=< 

H 

1 

Ds  =  D  +  ieD' 
D  =  {f -Fa)2-  k 
D'  =  f-Fa 


(31a) 


(31b) 


(31c) 


(31d) 


£  =  exp[  kz  +  i(ax  +  (3y)} 

£*  —  exp[fc(  -  i(a£  +  Prj)] 

In  the  special  case  F  =  0,  i.e.  for  wave  diffraction-radiation  without  forward  speed,  (31b)  yields 
A  =  1  if  k  =  /2,  i.e.  at  the  dispersion  curve  defined  by  the  dispersion  relation  D  =  0 .  The  flow 
related  to  the  Fourier  component  T  defined  by  (31)  can  be  analyzed  simply  and  effectively  using 
the  Fourier-Kochin  approach  expounded  in  Noblesse  and  Yang  [1] .  This  approach  is  used  below. 


Fourier-Kochin  representation  of  free-surface  effects 

By  substituting  (31a)  into  (16),  (21b),  (26b),  (29b),  (29c)  and  exchanging  the  order  in  which 
the  integrations  with  respect  to  the  Fourier  variables  (a ,  0)  and  the  space  variables  (x  ,y  ,z)  are 
performed,  we  can  express  the  free-surface  components  in  (21a)  and  (29a)  as 


47T  < 


kH 


un+u 


n 


n 


vn+v 


W7*  +  wj  j 


1  roc  roc 

=  lim  —  dp  da  ,  _ 
£ — *+0  7T  7— oo  7—oo  &  77, 


S4> 

f*£*  j  Su+kS j 
Sv+  k  Sj 
Sw+kS}u  J 


4-7T  < 


<f)H+<j)p  +  <t>r 

'  S4,+F2Sf  +  2rSf 

UH+  Up  +  u1/ 

1 

1  roc  roc  £* 

>  =  lim  -  d/3  da  —  < 
e— ++0  7 T  7—oo  7—oo  De 

SU+F2S£  +f2Sj 

VH  +Vp  +V^  +  V? 

Sv+F2Sp  +  f2Sj  +  2t  Sy 

WH  +  Wp  +  Wj?  +  wl? 

k  Sw+F2S$  +f2Sf  +  2 r  5“’ , 

(32a) 


(32b) 


if  F  =  0  or  F  >  0 ,  respectively.  The  Fourier  components  ,  uH  and  pp  ,  Up  in  the  repre¬ 
sentations  (26a)  and  (29a)  are  identical.  Indeed,  the  Fourier-Kochin  representation  of  steady  ship 
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waves  is  a  particular  case  of  the  Fourier-Kochin  representation  of  time-harmonic  ship  waves,  and 
is  obtained  by  setting  /  =  0  and  t  =  0  in  (32b) . 

The  spectrum  functions  S'*  Su,  Sv,  Sw  corresponding  to  the  components  <j)H,  uH  and  pH ,  uH 
defined  by  (16)  are  given  by 


s * 

Su 

>  =  f  dA  < 

sv 

J^HW 

Sw 

< 

\S 


(33a) 


u-n  —  (z  a  nx  + i  Pny+  knz)  <j) 

—i a  u  ■  n  —  i P  (ux  n)z  +  k  (ux n)y 
-i/3u-n  -  k  (uxn)x  +  i  a  {uxn)z 
[  ku  •  ft  +  ia  {uxn)y  —  i  ft  (ux n)x 

The  spectrum  functions  Sf ,  Sf ,  Sf  corresponding  to  the  components  uf  and  defined  by 
(21b)  and  (29c)  are  given  by 


w 

r 

\s> 

=  /  del 

Jthw 

Is/”  J 

< 

1  <t>  So 


(33b) 


ty 

-tx 

i(aty  —  Ptx)j  k 

with  So  =  exp[i(ax  +  f3y  )}.  The  spectrum  functions  Sf,  Sf ,  Sp ,  Sf  corresponding  to  the 
component  <f>p  ,  uf  defined  by  (26b)  are  given  by 


'sf 

—{ty)2  u  ■  u  +  txty  u  -t  —  iaty4> 

QU 

Op 

>  =  [  dd 

i  a  (ty)2  u-V  —  ia  txty  u  •  t 

> 

QV 

Op 

J  rHW 

ifi{ty)2u-  V  —  i(fitxty+  a)u  -t 

QW 

*F  J 

k  —  k  ( ty)2u •  V  +  ( ktxtv  +  a  P/k)u  -t  _ 

}So 


(33c) 


The  spectrum  functions  Sf ,  Sf ,  Sf  corresponding  to  the  component  ,  uf  defined  by  (29b) 
are  given  by 


f 

Sf 

ity4> 

Sf 

>  =  [  dC  < 

iu-t 

JpHW 

sf  ^ 

—(fi/k)  u  -t 

\ 

}So 


(33d) 


Similarly,  (31a)  and  (21c)  show  that  the  component  <j)p ,  up  in  (21a),  (26a),  (29a)  is  given  by 
the  Fourier-Kochin  representation 


4tt  < 


<t>p 

up 

vp 

wp 


1  roc  roc  c* 

)•  =  lim  —  d/3  da  — 
£—++0  7 T  7-00  7—00  D£ 


1 

— ia 

~ip 

k 


\Sp 


(34) 
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with  Sp  given  by  (35c) . 


The  next  section  shows  that  the  15  spectrum  functions  in  (32)  can  be  expressed  in  terms  of 
4  basic  spectrum  functions  S ,  Sf,  S* ,  SF.  These  4  basic  spectrum  functions  and  the  spectrum 
function  Sp  in  (34)  are  given  by 


S=[  dA[u  •  n  +  i%  (uxn)y  —  irr  (uxn)1]^ 
Jt.hw  k  k  J 


V 

i(f>(atv  —  f3tx)/k 

si 

>  =  f  dC< 

u  ■  t 

SF 

Jthw 

( txty+  a/3/k2)u-t  —  ( ty )2 u  ■  v 
<  * 

with  s=ekz  +  i(ax  +  0y)  £Q  =  +  0V) 


(35a) 


(35b) 


(35c) 

(35d) 


The  spectrum  function  S,  the  spectrum  functions  S?,  Sf ,  SF,  and  the  spectrum  function  Sp 
are  respectively  given  by  distributions  of  elementary  waves  over  the  boundary  surface  T,HW ,  the 
boundary  curve  YHW ,  and  the  mean  free  surface  EF.  The  spectrum  functions  S,  S* ,  SF  are 
defined  in  terms  of  the  normal  components  u-n,u-v  and  the  tangential  components  u  x  n ,  u-t  of 
the  velocity  u  at  the  boundary  surface  T,HW  and  the  boundary  curve  THW.  The  spectrum  function 
5*  is  defined  in  terms  of  the  potential  <f> . 


Transformations  of  spectrum  functions  in  Fourier-Kochin  representation 


Useful  transformations  of  the  spectrum  functions  in  (32)  are  given  below.  These  transformations 
are  performed  in  two  steps.  The  first  step  is  the  transformation  (37) ,  which  expresses  the  spectrum 
functions  S Su,  Sv,  Sw  given  by  (33a)  in  terms  of  the  three  spectrum  functions  S,  St,  S*  • 
The  spectrum  function  S ^  associated  with  the  potentials  4>n  and  (f>H  in  (32)  has  already  been 
considered  in  Noblesse  and  Yang  [1] .  Expressions  (33) ,  (29) ,  (36) ,  (34)  in  this  previous  study, 
where  the  vector  n  normal  to  the  boundary  surface  T,HW  points  outside  the  flow  domain  (instead 
of  inside  in  the  present  study) ,  yield  the  important  basic  transformation 


—  f  dA  (ianx  +  i  f3ny+  knz  )(j>  £  = 

JY.HW 

if  dA[r(uxn)y-r(uxn)x]£  +  SS 
Jzhw  lk  k 


with  Sf  given  by  (35b) .  The  spectrum  function  S  ^  in  (33a)  can  then  be  expressed  as 


S<P  =  S  +  S? 


(36a) 
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with  S  given  by  (35a) .  Expressions  (33a)  and  (35a)  immediately  yield 

< 


Sw=kS 


(36b) 


Section  A  of  Appendix  2  shows  that  (33a)  and  (35a)  also  yield 


Su+iaS  =  -ifiS*/k 


Sv+ifiS  =  iaS*/k 


(36c) 


Expressions  (36)  yield 


\ 

s + 

( 

1 

Su 

— ia 

►  =  i 

> 

sv 

-if3 

gw 

k 

\  / 

}S  +  S * 


(37a) 


with 


S*=< 


St 

-ifiSi/k 

iaSf/k 

0 


(37b) 


These  expressions  for  the  spectrum  functions  S*,  Su,  Sv,  Sw  —  with  S ,  St,  S*  given  by  (35)  — 
are  now  used  to  modify  the  spectrum  functions  in  the  Fourier-Kochin  representations  (32) . 

The  spectrum  functions  in  (32a) ,  corresponding  to  diffraction-radiation  of  time-harmonic  waves 
without  forward  speed,  can  be  expressed  as 


'  S* 

(  \ 
1 

su+ksy 

—ia 

>  =  < 

Sv+kS J 

-ip 

Sw+ kSf 

J  J 

k 

\ 

H  s  +  sf) 


(38a) 


with  S  and  St  given  by  (35) .  The  expression  for  the  spectrum  function  S $  follows  immediately 
from  (37) .  The  expression  for  the  spectrum  function  related  to  w  also  follows  immediately  from 
(37) ,  (33b)  and  (35b) .  The  expressions  for  the  spectrum  functions  related  to  u  and  v  are  verified 
in  Section  B  of  Appendix  2.  In  the  special  case  /  =  0,  corresponding  to  steady  ship  waves,  the 
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spectrum  functions  in  (32b)  can  be  expressed  as 


'  S^+F^S}  ' 

/ 

1 

Su+F2Sji 

>  —  < 

—ia 

SV+F2S£ 

-ip 

SW+F2S% 
r  / 

k 

\ 

(38b) 


Here,  S  and  SF  are  given  by  (35) ,  D  =  F2a2  —  k ,  and  5*  is  defined  by  (37b)  and  (35b) .  The 
expression  for  the  spectrum  function  related  to  w  follows  immediately  from  (37)  and  (33c) .  The 
expressions  for  the  spectrum  functions  related  to  4> ,  u  and  v  are  verified  in  Section  C  of  Appendix 
2.  In  the  general  case  fF  >  0 ,  corresponding  to  time-harmonic  ship  waves,  the  spectrum  functions 
in  (32b)  can  be  expressed  as 


'  S<t>+F2S$  +  2TSf 

\ 

1 

5u+F25>+/25V 

>  -  t 

—ia 

Sv+F2Svf  +f2S )  +  2  tSvt 

f  —  > 

-i(3 

k  Sw+F2Sf  +f2Sf  +  2 tS?  _ 

k 

\ 

where  D  =  ( Fa-f )2-  fc,  S*  is  given  by  (37b)  and  (35b) ,  and  S  is  defined  as 

S  =  S  +  £  S++  F2Sf- %  y  5* 
k  kk 


(38c) 


Expressions  (38c)  are  verified  in  Section  D  of  Appendix  2.  These  expressions  are  identical  to  (38b) 
in  the  special  case  /  =  0 . 


Fourier-Kochin  representation  of  waves 

Noblesse  and  Chen  [5,6]  and  Noblesse  et  al.  [7]  show  that  double  Fourier  integral  representations 
of  free-surface  effects  like  (32)  and  (34)  can  be  expressed  as  sums  of  wave  and  local  components: 


>  =  < 

'  <t>w+<t>L' 

u^+uy  —  up  J 

[uw+uL 

(39a) 


4>h  +  ^  +  <t>?  -  <j>p 

UH  +  Up  +  UF  +  UF 


} 

^  =  < 

1<PW+4'L\ 

ripj 

[uw+SL  ’ 

(39b) 


The  most  important  result  given  in  [5-7]  is  a  remarkably  simple  expression  for  the  wave  component 
(pw ,  uw .  This  expression  defines  the  wave  components  included  in  the  double  Fourier  integrals  (32) 
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and  (34)  in  terms  of  single  Fourier  integrals  along  the  dispersion  curves  defined  by  the  dispersion 
relation  D  =  0 .  It  follows  that  the  second  components  on  the  right  of  (38b)  and  (38c) ,  which 
involve  the  dispersion  function  D  as  a  factor,  do  not  contribute  to  the  waves  contained  in  the 


Fourier  representations  (32) .  Specifically,  the  wave  component  <j>w,  uw  is  given  by 


r  4>w ' 

( 

1 

47T  < 

uw 

VW 

>  =  -i  V)  fds 

D= 0 ^D=0 

sign  (£)')  +  ©  ^ 

—ia 

> 

-i(3 

w 

ww 

\  > 

k 

(40a) 


with  £*  —  eCfc-*(£a  +  f?/3)  as  given  by  (31d) .  The  identity  f2  =  k  at  a  dispersion  curve  D  —  0 ,  which 
follows  from  (31c) ,  was  used  in  (40a) .  (Da ,  Dp)  are  the  derivatives  of  the  dispersion  function  D 
with  respect  to  the  Fourier  variables  (a,/3).  D  and  D'  are  defined  by  (31c)  as  D  =  (f-Fa)2-k 
and  D'  =  f—Fa .  The  function  0  is  the  error  function 


e  =  erf 

where  a  «  1  is  a  positive  real  number.  The  wave-spectrum  function  Sw  in  (40a)  is  defined  by  (32) , 
(34) ,  (38)  as 


(40b) 


Thus,  waves  are  defined  in  terms  of  the  spectrum  functions  Sw  and  Sp  given  by  (40b)  and  (35) . 
In  the  special  case  /  =  0 ,  i.e.  for  steady  flow,  the  spectrum  function  Sw  is  given  in  terms  of  the 
spectrum  functions  S  and  SF,  which  do  not  involve  the  potential  (j>  and  are  defined  directly  in 
terms  of  the  velocity  u .  If  /  >  0 ,  Sw  involves  the  spectrum  function  Sf  which  is  defined  in  terms 
of  <fi ,  although  the  free-surface  boundary  condition  (30)  can  be  used  to  express  (p  in  terms  of  u . 


s  +  st 

- \ 

II 

o 

Sw=< 

S  +  F2Sf 

>-Sp  if  < 

/  =  o  ► 

'S+Zst+F’sr-lflS^ 

fF  >  0 
^  / 

The  representation  (40a)  of  the 


wave  component  (f>w ,  uw 

uw(o  =  ^w(i) 


agrees  with  the  relation 


used  in  Noblesse  et  al.  [8]  to  determine  the  velocity  uw  from  the  representation  of  the  potential 
<f)w  given  in  Noblesse  and  Yang  [1] .  Thus,  the  velocity  representation  obtained  in  the  present 
study  and  the  potential  representation  given  in  Noblesse  and  Yang  [1]  yield  identical  waves.  This 
agreement  needed  to  be  verified  since  the  relations  (38)  are  not  a-priori  obvious. 
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Conclusion 


A  new  fundamental  mathematical  representation  of  free-surface  flows  has  been  given.  This  flow 
representation,  called  velocity  representation,  is  given  by  (14)-(16)  for  generic  free-surface  flows 
associated  with  a  Green  function  G  of  the  type  defined  by  (8)  and  (12).  The  velocity  representation, 
and  the  classical  potential  representation,  are  given  by  (18),  (21),  (26),  (29)  with  (15)  and  (16)  for 
the  four  classes  of  flows  corresponding  to  the  particular  cases  of  (i)  an  infinitely  rigid  or  soft  free- 
surface  plane,  (ii)  diffraction-radiation  of  time-harmonic  water  waves  without  forward  speed,  (iii) 
steady  ship  waves,  and  (iv)  time-harmonic  ship  waves  (diffraction-radiation  with  forward  speed). 
The  potential  and  velocity  representations,  given  here  for  deep  water,  can  be  extended  to  uniform 
finite  water  depth,  and  indeed  can  be  extended  to  a  broader  class  of  problems  involving  dispersive 
waves  that  propagate  in  the  presence  of  a  planar  boundary. 

The  velocity  representation  defines  the  velocity  u  inside  a  flow  domain  in  terms  of  source  and 
vortex  distributions  with  strength  equal  to  the  normal  and  tangential  components  of  the  velocity  u 
at  the  boundary  surface.  Thus,  the  velocity  representation  does  not  involve  the  velocity  potential. 
This  property  is  an  important  difference  between  the  velocity  representation  and  the  classical 
potential  representation,  which  defines  the  velocity  potential  <j)  within  a  potential-flow  region  in 
terms  of  the  potential  (j)  and  its  normal  derivative  d<f)/dn  at  the  boundary  surface.  The  velocity 
representation  can  then  be  used  to  couple  a  nearfield  flow  calculation  method  based  on  the  Euler 
or  RANS  equations,  for  which  a  velocity  potential  cannot  be  defined,  and  a  farfield  potential  flow 
representation.  Another  notable  difference  between  the  potential  and  velocity  representations  is 
that  the  velocity  representation  defines  u  in  terms  of  first  derivatives  of  the  Green  function  G, 
whereas  u  can  only  be  obtained  from  the  potential  representation  via  either  numerical  or  analytical 
differentiation  of  </> .  Analytical  differentiation  of  the  potential  representation  involves  second-order 
derivatives  of  G ,  as  shown  in  (3c) . 

The  velocity  representation  provides  a  new  mathematical  basis  for  computing  free-surface  flows 
about  nonlifting  and  lifting  bodies  using  free-surface  Green  functions.  In  particular,  this  boundary- 
integral  representation  has  been  used  here  together  with  the  Fourier-Kochin  approach  expounded 
in  Noblesse  and  Yang  [1]  to  obtain  analytical  representations,  given  by  (40)  and  (35),  of  the  waves 
generated  by  an  arbitrary  boundary  velocity  distribution  for  the  three  basic  types  of  free-surface 
flows  defined  by  the  free-surface  boundary  conditions  (30),  i.e.  for  diffraction-radiation  of  water 
waves  without  forward  speed,  steady  ship  waves,  and  time-harmonic  ship  waves.  The  analytical 
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representations  of  waves  given  by  (40)  and  (35)  are  remarkably  simple.  The  Fourier-Kochin  rep¬ 
resentations  of  waves  obtained  here  from  the  velocity  representation  and  in  Noblesse  and  Yang  [1] 
from  the  potential  representation  are  in  agreement.  Practical  applications  of  the  Fourier-Kochin 
wave  representation  have  previously  been  reported,  for  the  case  of  steady  ship  waves,  in  Guillerm 
and  Alessandrini  [2]  and  Yang  et  al.  [3,4] . 

The  velocity  representation,  which  has  been  used  here  within  the  Fourier-Kochin  approach  to 
obtain  analytical  representations  of  waves  generated  by  a  boundary  velocity  distribution,  can  also 
be  used  to  obtain  corresponding  analytical  representations  of  nearfield  free-surface  flows.  Indeed, 
the  velocity  representation  given  in  the  present  study,  the  Fourier-Kochin  approach  expounded 
in  Noblesse  and  Yang  [1]  and  used  here,  and  the  practical  Fourier  representation  of  super  Green 
functions  given  in  Noblesse  and  Chen  [6]  yield  remarkably  simple  analytical  representations  of 
nearfield  flows  for  diffraction-radiation  of  water  waves  without  forward  speed,  steady  ship  waves, 
and  time-harmonic  ship  waves.  These  analytical  representations  of  nearfield  flows,  called  Rankine 
and  Fourier-Kochin  nearfield  flow  representations,  which  extend  the  Fourier-Kochin  wave  represen¬ 
tations  given  here,  will  be  given  elsewhere. 
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Appendix  1 

The  representation  (10c)  is  now  verified.  The  velocity  component  uF  in  (10a)  can  be  expressed 
in  the  form  (3a)  where  ux  is  defined  by  (3c)  and  (9)  as  /E  dA  V\  with 

( 4>  GF  +  cf>  Gfzz  )  nx  -  <(>  GFy  ny  -  cf>  GFZ  nz ' 

Vi  =  |  -4>GFX  nx  +  ( <j>GFz  +  <j>GFx)  ny  -<j>GFyz  n z  > 

^  <t>  Gfx  nx  +  (f>  GFy  nv  -  {4>GFX  +  4>GFy)nz . 

Expression  (10c)  for  uF  —  —  ux  holds  if  ux  can  be  expressed  as  /s  dA  V?  with 

{(j)x  ny  -  <j)y  nx)  Gy  -  ( 4>z  nx  -  4>x  nz)  GF 
V2=<  (d>y nz  —  <t>z ny) Gf  —  (4>x ny  —  <py nx) Gx  • 

~(<j>z  nx  -  (j)x  nz)  Gf  +  (<t>y  nz  -  (f>z  ny)  Gy  J 

i.e.  if  fj*  dA  (V \  —  V2)  is  null.  The  x,y,z  components  of  the  vector  Vi  —  V2  are  given  by 
Vx  •  n,Vy  ■  n,-Vz-  n  where  the  vectors  Vx,Vy,Vz  are  defined  by  (6).  The  divergence  theo¬ 
rem  applied  to  the  vector  fields  Vx,  Vy,Vz  and  the  closed  surface  E  completes  the  verification  of 
(10c) .  This  expression  is  at  variance  with  (19c)  in  Noblesse  et  al.  [8] ,  which  was  obtained  via  an 
incorrect  application  of  the  divergence  theorem  to  an  open  boundary  surface. 


Appendix  2 


A.  Expressions  (33a)  and  (35a)  define  the  function  Sf  in  (36c)  as 

Sf  =  J  ^wdA[ia  (uxn)x  +  i  /3  ( uxn)v+  k  ( uxn)z  ]  £ 

where  £  =  exp[ kz  +  i  (ax  +  fly)] .  We  have 

5*  =  [  dA  V£  •  (ux  n)  =  [  dA  n  ■  ( V£  x  u )  =  [  dA  n  ■  ( V  x  £u  ) 

Jxhw  K  1  Jxhw  v  y  Jehw  v  ' 

since  V  x  u  =  0 .  Stokes’  theorem  shows  that  5*  is  the  function  given  by  (35b) . 


B.  Expressions  (37)  yield 


r 

Su+kSf  =  -iaS -iflSl/k  +  kSf 
Sv+  kSJ  =  -i0S  +  iaSt/k  +  kSf 
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Thus,  one  needs  to  verify  the  relations 

Q  {-ipS*/k  +  kS?  +  iaS? 

\  iaS*/k  +  kSj+ipst 

where  So  =  exp[i (ax  +  fiy)].  Expressions  (33b)  and  (35b)  yield 

\~P 


-a 

II 

M 

> 

J  Jyhw 

Av 

<  j 

}S0 


L 


yHW 


dC{ 


AU 

Av 


a 


j^HwdC[u-t  +  i(af+ (3ty)(j)]£o 


We  have 


J^wdC[u-t  +  i(atx+pty  )</>]£ 0  =  j  dC  t-S7(</>£0)  =  0 

which  verifies  the  expressions  for  the  spectrum  functions  related  to  u  and  v  in  (38a) . 

C.  Expressions  (37)  yield 

S^+F2Sp  =  S  +  Sf  +F2Sp 
Su+F2Sp  =  —iaS  —  i/3  S*/  k  +F2Sp 
Sv+F2S£  =  -i(3S  +  iaS$/k  +F2S£ 

Thus,  one  needs  to  verify  the  relations 


0  = 


S$-SF+S?a2/k 


F2(S$-SF)  +  (l+D/k)S t 
F2(Sp  +iaSF )  —i  (1+D/  k)  St  P/  k 
F2(S£  +  ipSF)+i(l+D/k)Sta/k 

Here,  D  =  F2a2—  k .  These  relations  can  be  verified  using  (33c) ,  (35b) ,  and  (41) 

D.  Expressions  (38c)  and  (38b)  show  that  one  needs  to  verify  the  relations 


S%  +  iaSF-iSt{3a2/k2 
S$  +  ipSF+iSta3/k2 


2  tS$ 

1 

f2Sf 

—ia 

>  —  i 

f2Sf  +  2r  Sp 

-i(3 

f2sy  +  2tS™  j 

k 

\ 

f 

k 


2rp 

k  k 


2  ra  —f2 
k 


Thus,  one  needs  to  verify  the  relations 

kS*-aS?  +  pS*/k 
kSf  +  iaSt -ipSt/k 

k  S /  +  (2  F/f)  ks?  +  ips?  +  iast/k-i  (2  F/f)  k  S 
k  Sf  +  (2  F/f)  kSy-k  St  +  (2  F/f)  p  St 
These  relations  can  be  verified  using  (33b) ,  (33d) ,  (35b) ,  and  (41) . 


0  =  < 


(41) 
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